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Introduct ion 


The technique for determining the attitude of the magnetopause 
current layer which was introduced by Sonnerup and Cahill [1967] 
has been applied by several experimenters (Sonnerup and Cahill 
[1967, 1968], Cummings and Coleman [1968], Kauffman and Konradi 
[1969], Aubry et al. [1970]). However, the only surveys extending 
over large areas of the magnetopause have used data from Explorer 12, 
which provides measurements of the magnetic field components with a 
resolution of 24 gammas. The satellite spin improved the effective 
resolution in the spin plane; nevertheless, it appears desirable to 
further apply the technique to data with higher measurement resolution. 
This paper reports the study of magnetopause crossings of the Fifth 
Orbiting Geophysical Observatory (OGO-5), using data taken by the 
Goddard Space Flight Center magnetic field experiment which measures 
field components with a resolution of 1/4 gamma. The sampling rate 
of the three components is 1.7, 14, or 56 samples per second, depending 
on the mode of operation of the observatory. A total of 70 crossings 
were examined, of which 31 satisfied accuracy criteria of the calculation. 
These crossings were located at subsatellite local times between 1020 
and 0430, and solar magnet os pheric latitudes between -7° and 440°. 
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Applicatlon of the Technique 

The technique of Sonnerup and Cahill consists of determining 
the components of a unit vector n x which will minimize the variance 

N 

- {f I [Si • ft, - 1 • fij 3 

i=l 

where I?i represents one of N individual field measurements taken 

N 

U 1 r-» _ 

during a magnetopause crossing. The quantity B is equal to — ) Bj[ 

1 i-1 

n x is interpreted to be the normal to the magnetopause current 
layer. This computation, of minimization of the variance, is 
equivalent to that of finding the smallest principal axis of the 
variance ellipsoid defined by Equation [1]. The calculation also 
yields n a and n 3 , and <7 a and cr 3 , the directions and magnitudes of 
the other two principal axes. 

Each 0G0-5 crossing was analyzed as follows: The approximate 

time of crossing of the center of the current layer was determined 
by inspection. The measurements were then grouped in a series of 
intervals each centered on this time and having a successively 
increasing width (see Figure 1). The increment of data for each 
successive segment, and the largest Interval of data used were 
chosen for each crossing so that the innermost segments would 
contain only measurements taken inside the current layer, and the 
outermost would contain, in addition, magnet ospherlc and magnetosheath 
measurements. (By "inside the current layer", one means within a 
time interval when the field can be recognized as being transitional 
between the magnet ospherlc field and the magnetosheath field.) 


Typically, about one minute of data and ten or more nested segments 
were used. The three vectors n x , n a , and n 3 and their corresponding 
variances were determined for each segment. In practice, it was 
usually found that the direction of n x , n a , and n 3 for the inner 
segments varied considerably from one segment to another. As the 
segment width increased, a series of segments were commonly obtained 
which had relatively constant directions of n x . As still wider 
segments were used, the variance would increase and the direction of 
n x might again vary. The criterion adopted for an acceptable deter- 
mination of n x was that there be a consecutive set of segments, 
greater in number than those in the variable inner region, for which 
(1) the directions of n x all fell within a cone of total angle less 
than 15°; (2) the magnitude of was less than 0.3 times the magnitude 
of the vector change in the field across the magnetopause current layer; 
and, (3) a a /oj was <1/3. When these conditions were satlsified, 
a representative value of was chosen from one of the acceptable 
segment 8 . 

In fact, the average total cone angle for the 31 acceptable 
crossings was 6°. The data rate varied from 1.7 to 56 samples per 
second depending on the mode of operation of the observatory, but 
trial calculations showed that the calculation was insensitive to 
both the sampling rate and the exact time chosen for the center of 
the nested segments. The magnetopause crossings were found by 
scanning only the magnetic field experiment data; a relatively clear 
identification was required before an event was chosen for calculation. 
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The n x vectors were transformed into a spherical system, in 
which the polar axis is the earth-to-sun line and the zero direction 
fot the azimuthal angle is a line perpendicular to a plane containing 
the earth-sun line and the magnetic dipole axis of the earth (See 
Figure 2). 

A simple model surface, suggested by Mead (private communication) 
was used to generate theoretical normal vectors at the coordinates of 
the magnetopause crossings. The model has a geocentric radial distance 
proportional to sec 0/2 where 0 is the angle from the polar axis and is 
a surface of revolution about this axis. This surface flares out with 
increasing 0 slightly more than does the published Mead and Beard model 
[1964]. For example, the ratio of the radial distance of 0 = 90° to 
that at 0 = 0° is 4 % greater in this than in the published Mead and 
Beard model [op. cit,]. This general behavior is in accordance with 
experimental observations of magnetopause locations. 


Results 


(1) Deflection of the Current Laver 

Figure 3 is a plot of A0 versus 0 for oach acceptable 
calculation of . Here A0 Is the angle 0 of the experimentally 
determined n> vector minus the angle 0 of the n x vector obtained 
from the simple model. The abscissa, 0, is the 0 coordinate of 
the satellite position at the crossing. In general, A0 is seen to 
be greater for the magnet osheath -to -magnetosphere crossings than 
for the magnet osphere-to-magnetosheath crossings. 

Figure 3 is a plot of n^ versus 0, where n^ is the component 
of the experimentally determined unit vector n x perpendicular to 
the plane of constant 0 passing through the satellite position. 

The abscissa is the 0 coordinate of the satellite position at the 
crossing. The value of iy from the simple model would be zero 
everywhere, since the model is a surface of revolution about the 
polar axis. In this plot a positive value of n^ corresponds to 
the direction of increasing angle 0. The plot symbols have the 
same meaning as before. 

It will be noted that in this plot the two directions of 
crossing ere not clearly separated as they were in Figure 2. This 
tendency for the normal vector to be deflected to a larger or 
smaller value of 0, depending on whether the crossing is from the 
magnetosphere to the transition region or vice versa is consistent 
with a model in which the crossing occurs as the result of a 
propagation of a compression or expansion front over- the surface 
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of Che magnetosphere from the day to the night side. This model 
was suggested, for example , by Kauffman and Konradl (op. clt.) In 
their interpretation of Explorer 12 results. Aubry et al. (op. clt.), 
also Interpret their data to signify anti-solar propagating waves. 

If this model Is applicable, then it appears that a large majority 
of the present crossings must have occurred during the passage of 
such fronts past the satellite. 

The dotted lines in Figure 3 are least -squares fits to the two 
directions of crossing. The equations to the least square lines 
are: A9 * 17.7 - (0.21 + 0.15)9 for transition region to magneto- 

sphere crossings, and A9 ■ -19.9 - (0.09 + 0.17)0 for magnetosphere 
to transition region crossings. The units are degrees. The errors 
quoted for the slope are standard deviations, calculated under the 
assumption that a linear relationship exists between A0 and 9 for 
each type of crossing. These two lines converge slightly with 
increasing 0, but the convergence is not statistically significant. 
That is, to the accuracy of these data, there is no evidence of 
either an increase or a decrease of the deflection of the normal 
vector with increasing sun-earth-satellite angle, 9. 

(2) Magnitude of the Field Per p endicular to the Current Lave r 

In reconnection models of the magnetosphere, the magnitude of 
the magnetic field perpendicular to the magnetopause current layer 
is a measure of the reconnection rate, e.g. Levy et al. [1964]. In 
these models reconnection can only occur between the magnet ospheric 
and magnetosheath fields when they are inclined at an angle of 
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greater than 90° with respect to each other and then the perpendicular 
field will be Inward directed in the Northern Hemisphere and outward 
in the Southern Hemisphere. Levy et al. [op. cit.) predict that this 
perpendicular field will be 10-20% of the adjacent magnetospheric and 
magnetosheath fields when these latter fields are anti-parallel. When 
the magnetospheric and magnetosheath fields make an angle of less than 
90°, the perpendicular component is predicted to be small or zero. 

An attempt was made to test this model by using the 31 0G0-5 
calculations of magnetopause attitude. First, the angle between the 
transition region and magnet opsheric fields was determined by examining 
the data taken within one or two minutes on either side of the magneto- 
pause. The determination of this angle was, in many crossings, 
impossible on account of rapid direction changes in the transition 
region field. Of the 31 crossings, 17 were judged to be sufficiently 
stable to define whether the angle was greater or less than 90° (see 
Table 1). Ten were greater than 90°; they were all in the Northern 
Solar Magnetospheric Hemisphere. Seven were less than 90°; six of 
these were In the Northern Solar Magnetospheric Hemisphere, and one 
In the Southern Hemisphere. 


For each of the 17 crossings, the ratio of the normal component 
to the total magnetospheric field was computed, viz. (3 • n x ) /B^ 
where B M Is the magnitude of the magnetospheric field adjacent to the 
to the magnetopause and the average 3 * n x was taken over the segment 
of data used to determine n x . A standard deviation cr x /B^ was also 
computed for each ratio where o x , defined In Equation 1, was obtained 




•>-v. 
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from the same segment of data. The average of the ratios (If * n 1 )/B| k | 


for the ten crossings with fields at greater than 90 wcs -0.03 4 +.09, 

end for the seven cases where the fields were at less than 90°, the 

ratio was -0.05a +.06. Here, a negative sign corresponds to an inward 

directed field. The errors quoted are standard deviations computed 

from the set of ten (seven) values of the ratio R ■ (If * ft, )/Bvi using 

N * 


the formula: standard deviation 


the value of R for the j** 1 crossing 


-V[(R4 - R) a /(N - 1)]* where is 

A- ^ 

R - 2/ R j/ N > 


and N - 10 (or 7). 


As shown in Table 1 the standard deviation c^/B^ of the quantity 
Rj for any single crossing was generally smaller than the standard 
deviations computed for each of the two averaged values of R; that is, 
there were statistically significant variations in Rj from one crossing 
to another; individual crossings yielded a range of values of Rj from 
approximately +0.1 to -0.2. Because of this variation of Rj from one 
crossing to another, and since both positive and negative components of 
the magnetic field normal to the magnetopause surface were observed when 
the adjacent magnet ospherlc and magnetosheath fields are anti-parallel , 
the average ratio, R, taken over a number of crossings is not a good 
parameter for determining reconnection rates in this area of the 


magnetopause . 
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Concluslons 

(1) The average ratio of the normal component of the magnetopause 
field to the adjacent magnet oapheric field was 3+9 percent for 
ten 0G0-5 crossings which had occurred when the adjacent magneto - 
spheric and transition region fields were inclined to each other 
at an angle of greater than 90°. The average rat^c was 5+6 
percent for seven crossings when the fields were inclined at an 
angle of less than 90°. The errors quoted are largely the result 
of the variability of the ratio from one crossing to another. This 
variability precludes, in this region of the magnetosphere, the use 
of the average value to accurately test reconnection models that 
predict ratios of ten percent or less. 

(2) Straight line fits to the plots of deflection of the magnetopause 
current layer versus the sun -earth -satellite angle for each direction 
of crossing show no evidence of a growth of decay of the deflection as 
this angle Increases. 

(3) The direction of these deflections is consistent with a model in 
which the 0G0-5 crossings occvr during the passage of compression or 
expansion fronts propagating from the day to the night side of the 
magnetopause. This same interpretation was first proposed by Kauffman 
and Konradi (op. clt.) for Explorer 12 magnetopause crossings. 


ACKNOWLEDGMENTS 


The assistance of H. J. Gillis, who performed the mathematical 
programming, is gratefully acknowledged. 


REFERENCES 


Aubry, Michel P„, Margaret G. Kivelson, and Christopher T. Russell, 
Motion and Structure of the Magnetopause, Publication No. 865 . 
Institute of Geophysics and Planetary Physics . U.C.L.A., 1970. 

Cummings, W. D . , and R. J. Coleman, Jr., Magnetic fields in the 
magnetopause and vicinity at synchronous altitude, J. Gepphys. 

Res . . 73, 5699-5718, 1968. 

Kauffman, Richard L., and Andrei Konradi, Explorer 12 magnetopause 
observations; large scale non-uniform motion, J. Geophys. Res. . 

74, 3609-3627, 1969. 

Levy, R. H. , H. E. Petschek, and G, L. Siscoe, Aerodynamic Aspects 
of the Magnet ospheric Flow, AIAA Journal . 2, 2065-2076, 1964. 

Mead, Gilbert D., and David B. Beard, Shape of the geomagnetic field 
solar wind boundary, J. Geophys. Res. . 69., 1169-1179, 1964. 

Sonnerup, B.V.O., and L. J. Cahill, Jr., Magnetopause structure 

and attitude from Explorer 12 observations, J. Geophys » Res.. 72, 
171-183, 1967. 

Sonnerup, B.V.O., and L. J. Cahill, Jr., Explorer 12 observations 

of the magnetopause current layer, J. Geophys. Res. . 73 . 1757-1770, 
1968. 



FIGURE CAPTIONS 


FIGURE 1. 


FIGURE 2. 


FIGURE 3, 


FIGURE 4. 


N segments of successively increasing width, centered 
on the magnetopause current layer. 

The geocentric spherical coordinate system used. 

The polar axis points to the sun; the zero direction 
of the azimuthal angle is a line perpendicular both 
to this polar axis and to the geomagnetic dipole axis. 

A plot of A0 versus 9. A9 is the difference between 
the 9 coordinate of the experimentally determined n x 
vector and the 9 coordinate obtained from the simple 
model surface. The abscissa 9 is the 9 coordinate of 
the satellite's position. The symbol X is used for 
magnetosheath to magnetosphere crossings and the symbol 
• for magnetosphere to magnetosheath crossings. 

A plot of n^ , the 0 component of the experimentally 
determined n x vector, versus the 0 coordinate of the 
satellite's position. Note that the simple model 
surface would predict that n^ would be zero everywhere. 
The plot symbols have the same meaning as in Figure 3. 



POLAR AXIS 



X MAGNETOSHEATH TO MAGNETOSPHERE 
• MAGNETOSPHERE TO MAGNETOSHEATH 



FIG 3 


A COMPONENT OF 


X MAGNETOSHEATH TO MAGNETOSPHERE 
• MAGNETOSPHERE TO MAGNETOSHEATH 
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FIG 4 




TABLES 


TABLE 1. 


The 31 acceptable calculations are summarized here. 

A -fve sign in the column (B • n 1 )/B^ denotes an 
outward directed normal component. The column headed 
"subtended angle" refers to the angle between the 
magnetospheric and magnetosheath fields adjacent to 
the magnetopause. 
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TABLE 1 : Summary of Calculati 
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